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Abstract 

We show that the generating functional for hard thermal loops with external gluons in QCD 
is essentially given by the eikonal for a Chern-Simons gauge theory. This action, determined 
essentially by gauge invariance arguments, also gives an efficient way of obtaining the hard 
thermal loop contributions without the more involved calculation of Feynman diagrams. 
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The Chern-Simons (C-S) action made its appearance in physics hterature over ten 
years ago as a mass term for gauge fields in three dimensions ^. Studies since then have 
revealed many interesting properties of this action. The Abelian version can be used 
for spin transmutation, converting spin zero bosons into anyons, for example ^. The 
correlators of Wilson lines in a pure C-S theory are related to the the polynomial invariants 
of knot theory ^. Pure C-S theory is also closely related to conformal field theory and 
the Wess-Zumino-Witten (WZW) action in two dimensions Also actions related to 
the C-S action can be used for self-dual gauge fields and integrable systems ^. Finally 
there is an intriguing class of vortex solutions in spontaneously broken C-S theory ^. 
However, despite this bounty of interesting results there have not been many realistic 
physical systems for which the C-S action is relevant. In this letter we show that the 
C-S action, more precisely its eikonal, is part of the effective action for describing the 
gluon plasma in Quantum Chromodynamics (QCD). This action, determined essentially 
by gauge invariance requirements, gives an efficient way of obtaining the hard thermal loop 
contributions, without having to calculate the corresponding Feynman diagrams. The C-S 
connection is particularly interesting in view of the possibility of producing the quark-gluon 
plasma in heavy ion collisions in the near future. 

We consider QCD at temperatures well into the deconfinement phase; i.e. we have a 
'hot' plasma of gluons. The effective action mentioned above is more precisely defined as 
follows. Pisarski has shown that a partial resummation of Feynman diagrams in thermal 
QCD is necessary to obtain gauge invariant results, for example for the gluon damping 
rate in the plasma ^. The resummation amounts to the following. We calculate the one- 
loop diagrams of thermal QCD; the relevant kinematical regime corresponds to the loop 
momentum being much larger than the external momenta. These are the so called hard 
thermal loop contributions. For these, the external momenta are typically of the order of 
gT where g is the coupling constant and T is the temperature; the loop momentum being 
hard, i.e. at least of the order of T, is the region of interest. The leading contributions 
are proportional to T^. The generating functional for hard thermal loops is the effective 
action. Thus, once the high temperature contributions of the hard thermal loops have 
been obtained, calculations can be done for any process starting from the effective action. 
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The result is gauge invariant and consistently accounts for all terms of a given order in 
coupling constant ^. 

Many authors ^'^ have written down versions of this generating functional. The results 
involve a null vector = (1,Q) and integration J dfl over the directions of the unit 
vector Q. Diagrammatically, Q arises as follows. Thermal loops describe the absorption 
and emission of particles from the surrounding medium or thermal bath. These particles 
are on mass-shell and thus the loop integration is only integration over the momentum 
three-vector with a distribution of the Bose-Einstein or Fermi-Dirac form. Integration 
over the magnitude of this vector is then carried out leaving the angles, described by 
the unit vector Q, for the final integration. We shall also need the following coordinates 
{u, V, xt), u = |Q' ■ X, V = ■ X, Q ■ xt = 0, where = (1, —Q). The components of 
the gauge field along Q wiU be denoted by A"-, viz. = -it°'A'^, Q ■ A°- = A"-, Q ■ A = A. 
{t"} are a basis of the Lie algebra of the group G, chosen as hermitian matrices in the 
fundamental representation with Tr(t"t^) = ^S""^. Following reference 8, the generating 
functional for hard thermal loops with external gluons has the structure 

Utt 

C = Cq for quark loop contributions and C = Cq for gluon-ghost loop contributions; Cq 
and Cg are the quadratic Casimirs for the quark and adjoint representations respectively. 
We shall not display the coupling constant g in what follows, as it can be recovered by 
gAf^ and an overall The first term in brackets is the well known mass term for 
the time-component of the gauge field. The second term may be considered as what is 
necessary to render T gauge invariant. The information given by the diagrammatic analysis 
of hard thermal loops is that it can be written as J dfl W, where W is a functional of 

The condition for the gauge invar iance of F in (1) is 

j dnSW = A'K j dt^x Alu" (2) 
where SA^^ = d^uj + [^^,a;]. Equation (2) is realized by 

5W ^ I d'^x A'^uj" (3) 



■(27r j d^x A^qA^q + J dnW) (1) 



One can check that (3) is indeed the way gauge invariance is reahzed, again by analysis 
of diagrams. Thus it is equation (3) that we must solve. We rewrite this, using the 
transformation law for A"', as 



where / = ^ + ^A. Our analysis so far parallels reference 8. However, we shall now solve 
(3,4) in terms of the eikonal for the Chern-Simons theory. 

We begin by briefly recalling some aspects of the pure Chern-Simons theory Con- 
sider the C-S action 

S=^ Jd^'x Tria^d^aa + |a^a,aa)e'^"" (5) 

Here is the Lie algebra valued gauge potential, = —it°'a'^^. We shall use complex 
coordinates z^z^ z = x + iy, to describe the spatial dimensions. The time-component oq 
can be set to zero as a gauge choice. The equations of motion then tell us that a^, az are 
independent of time but satisfy the constraint 

dzas - dsa^ + [a^, a^] = 0. (6) 

This can be solved for function of a^, at least as a power series in a^. The result is 



_\^(_-t\n-l f ^ ^1 d Ctzjzi, Zi) . . .azjZn, Zn) 



This can be checked easily using dz^^ — 7t5^'^\z — z'). 
Define the functional /[a^] such that 

51=!^. J d^xTr{az[az] Saz) (8) 

Since is conjugate to a^, / so defined is the action evaluated for the classical motion; thus 
it is Hamilton's principal function or the eikonal for the C-S action. (Recall that the eikonal 
for one-dimensional particle mechanics is the integral of pdx, where p is the canonical 
momentum, specified as a function of x by fixing the energy. We have an analogous 
situation where is specified in terms of by equation (6).) We can write / as 

J _.j^^Y^ (-1)" f d^zi d'^Zn Tr{az{zi,zi) . . .az{zn,Zn)) 

^ n J TT TT Zi2Z23---Znl 
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where Zij = zi — Zj. (This result is fairly standard; for an elegant way of obtaining and 
rewriting this, see reference 9.) / is in fact the WZW action. If we parametrize the two- 
dimensional gauge field as = —dzU where U is valued in the complexification of 
the group G, (9) can be written in the more conventional form of the WZW action as 
/ = —ikSwzWi where 

M3 

(10) 

(As usual = X [0,1] with U{z,z,0) = 1, U{z,z,l) = U{z,z).) The WZW action 
is thus the eikonal for the C-S action. For our discussion below, will not be just a 
two-dimensional gauge field. Equation (9) will thus be the more useful form. Also, since 
k is not relevant for our discussion, we shall henceforth set it to one. 

Returning to equation(4), notice that it is of the form of a zero curvature condition. 
We first do a Wick rotation to Euclidian space so that 2u z, 2f — > z, du '2.8 z-, 
dy 2dz. Defining = — / and = ^A, equation (4) is seen to be identical to (6). 
Hence the solution for W is given by 

W = -^Jd'^x A" A" - AmI[A/2] (11) 

where, from (9) and the definition of a^. 

The potentials now depend on all four coordinates; however since (4) does not involve 
differentiations with respect to the transverse coordinates xt-, all potentials in (12) have 
the same argument for these coordinates. In other words, the transverse coordinates in 
(12) only play the role of parameters on which the ^'s depend. Using (11) in (1) we have 
the generating functional in terms of the eikonal /, 
„ CT2 



127r 



( j d^x {2tt A'^A^ -J J dn - 47rz J dfl I[A/2] 



(13) 



We have not actually calculated Feynman diagrams to arrive at (13). The only in- 
put from a diagrammatic analysis has been the structure J dQ W. However it is easy to 



check that the n-point functions calculated from (13) agree with the explicit diagrammatic 
evaluation of hard thermal loops. The n-point functions in momentum space are given by 



The two-point function is given by 



5Al\{x^)...5All{xn) 



(14) 

A=0 



r«t(xi,X2) = 5'^'^^ [47r5^o5.o5(^n^i - ^2) 

-I j dn Q,Q.{S^^\x^ - X,) - ^^^^fi-^)] (15) 

We need the Fourier transform to obtain the expression in momentum space. This is 
straightforward for the first two terms. For the last term we have an expression of the 
form 

f d'^Z2 h{z2,Z2) „ f d'^Z2 h{z2,Z2) 

H= / 73 ^ = -5^1 / — — (16) 

J -K {Zi- Z2Y J TT (21 - Z2) 

where h is an exponential of the form exjpi{kzZ + kzz). Using d^j^^-^ = 'k5^'^\z — z'), we 
get 

d,,H = -d,,h (17) 

This can be easily solved for H. The Fourier transform of (15), after Wick rotation to 
Minkowski space, with 2kg k • Q, 2k z k ■ Q', gives 

K'^^ = S'''j^{^^S^oS.o-M (18) 

where ff,^, = J dfi Q^Q,y^. 

The three-point function involves the factor (^12^23^31)"^ in addition to the transverse 
5-function and color and factors. Using the splitting 

' (^-^) (19) 



^12223^31 (^12)^ ^13 ^23 

the Fourier transform can be evaluated by the same method as in (16,17) to obtain 

r*. = f^/d. Q,Q.Q. - ^) (20) 
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Expressions (18,20) agree with the diagrammatic evaluation of hard thermal loops 
We have checked the four-point function in a similar way. For this and for the higher point 
functions, a splitting formula analogous to (19) is very useful. It is given by the Ward 
identity or in other words, by the recursive buildup of the correlator of currents in the 
WZW model. Combined with the Fourier transform method we have used, this gives an 
efficient way of calculating the higher point functions. 

The vectors Q, Q', in terms of the coordinates u, v, define a two-dimensional subspace 
in spacetime. Our results indicate that at high temperature, as far as the hard thermal 
loops are concerned, the dynamics is essentially the C-S dynamics for the components 
AuiAy. The final results, of course, do not depend on the choice of this subspace since we 

— * 

integrate over the orientations of Q. The choice of this subspace and the integration over 
the orientations of Q can be incorporated into the action by using, instead of (5), 



Q ~ Q'-plane. This action is similar to the Kahler- Cher n- Simons (KCS) action considered 
in reference 5. The difference is that for us u, being restricted to directions transverse to 
the Q — Q'-plane, is degenerate. As for the KCS theory, the equations of motion tell us 
that the fields do not depend on the extra fifth dimension in the action (21). 

Since the C-S action is odd under parity, its presence in a QCD calculation may be 
potentially worrisome. However, we do not have any parity violation because all our results 
are integrated over the orientations of Q. Only the parity preserving contributions to the 
n-point functions survive this integration. Also for the non-Abelian C-S action, one has to 
address the issue of quantization of the coefficient of the action. Again, the integration over 

— * 

the orientations of Q shows that there is no quantization. The quantization arises, in the 
usual analysis, from the requirement of invariance under homotopically nontrivial gauge 
transformations. In our case, there are no nontrivial gauge transformations consistent 
with the angular symmetry imposed by the integration over the orientations of Q. In 
other words, the relevant winding number corresponding to maps of the three-sphere into 
the group cannot be defined in a way that is invariant under the Q-integration. 




(21) 



where a;, 



|(m^nj^ — m^rii^), with m^,nj^ defining a basis for vectors transverse to the 
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We conclude by rewriting the effective action in anotlier way. Defining = |Q' 
A, ttz = \Q ■ A, wliere is no longer related to by equation (6), we can write 

d^Q 1 



K[a,,a,] = -2C 



^ j d^x Tr(az%) + il + il 



(22) 



Here N{Q) is the Bose-Einstein distribution and / is given by / (from (12)) with the 
change — > a^, z ^ z. The density K is gauge invariant; apart from integration over 
the tranverse coordinates, it can be written as Tr \og{DzDz) where D^, Dg are the cor- 
responding two-dimensional covariant derivatives in the adjoint representation.^^ It can 
also be considered as the Kahler potential associated to the symplectic structure for the 
C-S action. ^'^ Following reference 13, it may also be possible to relate this to the forward 
scattering amplitude for high energy gluons on a gauge field background. These issues will 
be discussed in more detail elsewhere. 



We thank A.Mueller, J.Schiff and E.Weinberg for discussions. 
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